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ATTAINABLE CIRCULATION ABOUT AIRFOILS IN CASCADE 


By Arthur W. Goldstedt and Artur Mager 


SUMMARY 

From, consideration of available information on boundary- 
layer behavior, a relation among profile thickness, maximum 
surface velocity, Reynolds number, velocity diagram, and 
solidity is established for a cascade of airfoils immersed in a 
two-dimensional incompressible fluid flow. Several cascades 
are computed to show the effect of various cascade design para- 
meters on minimum required cascade solidity. Comparisons 
with experimentally determined blade performance show that 
the derived blade loadings are equal or higher for moderate flow 
deceleration and somewhat lower for large deceleration. Blades 
with completely laminar flow appear practical for impulse or 
reaction blading. 

INTRODUCTION" 

In the design of compressor and turbine blading, the choice 
of blade sections and spacing is usually made on the basis of 
experimental results from several cascades. An enormous 
amount of data is required to cover the entire range of inlet 
and outlet angles, Reynolds numbers, blade sections, spac- 
ings, and maximum surface velocities. 

The relations among blade circulation and blade thickness, 
Reynolds number, maximum surface velocity, and spacing 
were determined at the NACA Lewis laboratory during 1949 
in a qualitative fashion for two-dimensional, incompressible 
flow. This problem has been studied empirically (references 1 
and 2) without any attempt to relate the problem to the 
basic determining factor — the boundary layer. The analysis 
presented herein establishes such a relation on a more rational 
basis. 

As the gas flows through a cascade of airfoils, the inlet 
velocity on the suction surface of a blade is raised to some 
maximum local value and then decreased until, near the 
trailing edge, its value is in the neighborhood of the exit 
velocity. A blade cascade will operate in a compressible 
fluid with less likelihood of local shocks, resultant losses, and 
choking if the design is such that the maximum surface 
velocity is maintained as low as possible. The preassigned, 
low maximum surface velocity should extend over as large a 
portion of the suction surface as possible before diffusion in 
order to raise the circulation to a maximum. Because the 
velocity must decrease near the trailing edge from its maxi- 
mum value to a value close to the exit velocity, the maximum 
possible rate of diffusion will, in effect, determine the arc 
length on which the suction-surface velocity is kept at its 


maximum. The blade circulation attainable is therefore 
determined by the maximum diffusion rate possible. 

In dealing with incompressible flow, the same general 
type of velocity distribution is specified herein, and the 
diffusion rate is determined by several conditions to avoid 
separation of flow from the blade surface. Once the possible 
velocity distribution on the suction surface is determined, 
an approximate computation is used with a specified velocity 
diagram and blade thickness to compute the pressure- 
surface-velocity distribution, blade circulation, and spacing. 

The results reported herein are not to be interpreted as 
giving an exact velocity distribution of airfoil sections for 
application. In the first place, the assumption of no separa- 
tion in some cases reduces obtainable lift and increases the 
blade solidity and the surface for friction, which increases 
skin friction loss. If a small region of separated flow is 
permitted, the increase in form drag may be insufficient to 
increase total drag in view of the increase in blade lift and 
decrease in blade number and skin friction. Another limita- 
tion is possible in the applicability of the boundary-layer 
equations used in prescribing the suction-surface-velocity 
distribution. The velocities deduced from these equations 
are unconventional in that the form factor of the boundary 
layer (ratio of displacement to momentum thickness) 
rapidly rises to the maximum permissible value and then 
remains nearly constant; whereas with conventional velpcity 
distributions, it slowly increases and then very rapidly rises 
near the separation region. Thus the empirical equations 
describing the development of the boundary-layer form factor 
may not be accurately applicable for such velocity distribu- 
tions. In order to demonstrate the several solutions that 
might be applicable, the surface velocities were prescribed 
according to three rules. Gamer’s equations (reference 3) 
were used with the form factor first rising rapidly and then 
remaining nearly constant with limiting values of either 1.84 
or 2.14. Also a diffusion rate was chosen equal to 0.885 of 
the value prescribed by Kalikhman (reference 4) for 
separation. 

A final consideration lies in the approximate nature of the 
procedure for computing the lower surface-velocity distribu- 
tion. In this process, the effect of thickness is assumed to 
result in a surface-velocity component, which varies as a 
second-degree parabolic function of the surface arc length. 
The results consequently apply to a special family of thick- 
ness distributions, which may not be the optimum for the 
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prescribed suction-surface-velocity distribution. Moreover, 
the estimate of the effect of thickness is inexact so that the 
estimate of the airfoil-thickness distribution will differ 
somewhat from the actual thickness of an airfoil with the 
prescribed surface velocities. _ 

The absolute values of blade circulations are inexact; some 
very definite trends are nevertheless demonstrated, particu- 
larly with respect to maximum surface velocity, Reynolds 
number, and blade thickness. 

SELECTION OF VELOCITY DISTRIBUTION ON SUCTION 
SURFACE OF BLADE 

GENERAL CONSIDERATIONS 

As mentioned in the “Introduction, ’’ the velocity distribu- 
tion on the suction surface to be investigated will consist of a 
velocity that is supposed to rise to its maximum value in a 
negligibly short distance from the leading-edge stagnation 
point and to maintain that value as far back along the airfoil 
as possible. The velocity must then decrease from the 
maximum value to the value near the trailing edge that differs 
slightly from the downstream velocity. A deceleration with 
a maximum safe velocity gradient without flow separation 
between the maximum velocity and a final velocity is re- 
quired. In order to find the suction-surface- velocity dis- 
tribution, it is necessary to investigate the development of 
the boundary layer, which determines the allowable rate of 
diffusion. The general shape of the velocity-distribution 
curves is shown in figure 1 , which also indicates some of the 
nomenclature, (All symbols used are defined in appendix A.) 



LAMINAR BOUNDARY LAYER 

Region of constant velocity. — For the region of constant 
velocity U(x) = U max , the equation for the momentum thick- 
ness of the laminar boundary layer is 

0=0.664 (1) 

V u max 

where 

0 momentum thickness of boundary layer 
v kinematic viscosity of gas 
x arc length on airfoil suction surface 
U suction-surface velocity just outside boundary layer 
The subscript max indicates maximum value of surface 


velocity. The Reynolds number of the boundary layer 
based on momentum thickness is defined as 


R^=Ud[v ( 2 ) 

The form factor II is defined as 

Hs5*/0 • T3) 


where 5* is the displacement thickness of the boundary 
layer. In the region of constant velocity, 17=2.614. 

Region of diffusion. — For diffusion with a laminar bound- 
ary layer, the approximate method of Loitsianskii (reference 5) 
was used. According to that system of equations, sepa- 
ration will occur if the proportional change in velocity dUjU 
per unit of arc length dx equal to the generalized momentum 
thickness L reaches a value of —0.08884. That is 


LdU 
U dx 


y, lj au 


where 


U 


- L^eR=e 2 - 

V 


( 4 ) 

( 5 ) 


must be greater than —0.08884. A value of D L =~ 0.06618 
was arbitrarily selected to give a high diffusion rate without 
introducing laminar separation. Loitsianskii’s equation for 
boundary-layer growth then becomes 

D l fd?U\ 

1 dU 

ulx^~Jdtrj 

where F, a function of D L , is a constant equal to 0.8517 for 
D l =- 0.06618. 

This equation is integrated to give 


^f-—KU~ FtVL 

dx 

By substituting for dU]dx from the definition of Z? t , the con- 
stant of integration K is evaluated to give 



Dl 

51+4- 


<£T 


— 0.08425 


(6a) 


where the subscript l indicates the value at the beginning of 
the laminar diffusion process. When equation (6a) is 
differentiated and 7 (7/ 17 is eliminated by means of equation (4), 
there is obtained 


, dx dU D l dL 


which is. integrated to 


L 


-=1+(D l +F) (jxf)= 1 + 0.7855 


(6b) 


Equations (1) and (5) are used to relate the constants x t 
and Li 

Li=0.441xj 
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From equation (5) 


Z £ fi /9YZ7 
Li di R, \6i) U i 


or 


also 


0__(L_ ZWFEJ _/i Y’ 5421 

er\L' u) -\Lj ~\L t ) 

F+2Dl 

12/ L\2if+Dl) _/ Z\°‘ 4ii79 

srVrJ 


The velocity potential defined by 

Udx 


is found from 


<p— •P i 

L t U, 


,-r 

Jo 

j ~rir;JL [(zD ' i_1 ] 


(6c) 


(6d) 


f6e) 


The form factor of the boundary layer H in the region 
D l = — 0.06618 has a value of 3.214 (reference 5). 

In applications, it is convenient to know the position x t in 
percentage of the total length of the suction surface A'. 
From equation (6b), 


AT_££_ 1 

Zf-Z+F+2?, 



(7) 


The subscript/ applies to the final value at the trailing edge 

Jj. 

of the blade. In terms of velocities, with y) = 044 l’ 


£=2.26S+ ? _ 


■D] 


f+d l 


Then 


% 1 ' 1 


X 

1 

, , 0.441 / 

U r \ d l 

1 1 f+d l \ 

u-J _1 J 


(7a) 

( 8 ) 


The constant <p ifUjjLi is computed from the equations for 
the nondiffusing section giving 

j^-2.268 (61) 


In figure 2, the variation of U/Ui and R/R t with (x—Xi)/L t 
is shown for diffusion with the laminar boundary layer. In 
applying these curves for any airfoil, the scale of the diffusing 
region must be adjusted for the value of L t that is equal to 



Figusj: 2.— Development of velocity and laminar boundary layer with constant diffusion coefficient. Z>l= -0.06618; 12=3.214. 
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Velocity ratio, I/jr/U*,* 

Figure 3.— Relation between sactton-aurfase potential and velocity ratio for laminar diffusion 
with Di.~-0.OGM8. 


0.441 x t . In figure 3, the curve, for the entire suction-surface 
potential is shown with X as the unit of length rather than 

u 

TRANSITION 

As the air flows over the blade surface, the thickness of the v 
laminar boundary layer will increase until a point is reached 
where the character of the boundary-layer flow will change 
from laminar to turbulent. The boundary layer will then 
develop according to different laws and will, therefore also 
sustain a different pressure gradient without separation. 

Transition occurs at a certain value of R tr , which will 
depend on the free-stream turbulence intensity and scale, the 
pressure gradient, the roughness and the curvature of . the 
surface, and on heat transfer through the boundary layer.. 
An R tr value of 2.50 has beim selected, .for this analysis in 
consideration of the high level of free-stream turbulence in 
turbomachinery, which tends to reduce R* r to a low value. 
According to GruschwiU (reference 6), observed boundary 
layers were always laminar if 2?<l250 regardless of the 
intensity of free-stream turbulence, and turbulent boundary 
layers have been observed for ff>250. Gruschwitz neglects 
to mention the pressure gradient at which these observations 
were made and as a consequence the assumption j? (r =250 
might be somewhat low for flow's with zero pressure gradient. 
It should be noted that if for any reason transition is delayed 
so that the rapid diffusion rate prescribed for turbulent 
boundary layer is applied to a laminar layer, separation 
results. Conversely, the effect of pressure gradient is to 
reduce markedly the value of A ir ; even though delayed in 
the zero pressure-gradient region, the transition will there- 
fore probably occur at the inception of the diffusion process. 


TURBULENT BOUNDARY LAYER 


The momentum equation for the boundary layer is 


6 dU 
dx 


: + T7W (2+iI) -pir 8 


where 

t surface shear stress 
p gas density 

If the definitions for friction coefficient 


( 9 ) 



GO) 


generalized momentum thickness 
'V . . .. - T=8R' /s 




( 11 ) 


and diffusion coefficient 



D r = tt 


TdU 


U dx 


( 12 ) 


are substituted into the momentum equation as in reference 3, 
there results 



dT = l 

dx 6 


[/-*(*+?)] 


(13) 


In order to complete the solutions of this equation for arbi- 
trary values of D T , two more equations are required: one to 
determine values of/, and one to determine the development 
of. the form factor II. 

Region of constant velocity. — In the region of constant 
velocity D T = 0, so that the momentum equation reduces to 


dT 7 , 
dx 6 J 


(14) 


Ealkner’s data (reference 7) for surface shear w’ith constant 
velocity show that 


*" /= 0-006534 (15) 

Equation (14) may therefore be integrated to give 


/(*-*„) (16) 

where the subscript tr indicates the value at transition. 

Conditions for diffusion without separation; — The system 
of equations proposed by Garner (reference 3) was selected 
to describe the development of the turbulent boundary layer 
because of its simplicity. In selecting a law for friction, 
Gamer _examincs data for 7?>1250 and shows that / is 
constant with a value of 0.006534 if — 0.01<C lb <(0.004. 
(This is. Falkner’s value of / for D r =0, reference 7.) No 
correlation is obtained for D r < — 0.01 . In consideration of 
the lack of correlated data for high diffusion rates, 0>D T > 
— 0.01 was selected for the range of diffusion rate and there- 
fore /= 0.006534. is used for the friction coefficient. 
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The development of the form factor H required for 
solution of the momentum equation is described by an em- 
pirical equation developed by Garner in reference 3. 
The relation is 

T*!§= (17) 
where .4=0.0135. 

The problem is then to define a velocity distribution that 
will avoid separation of the turbulent boundary layer. Ex- 
perience shows that when H reaches a value of about 2.0 
with conventional velocity distributions it subsequently rises 
rapidly and separation occurs. Values of H for separation 
are therefore usually given in the range 2.0<//<2.6. 
It is therefore proposed that one possible criterion for 
prescribing the velocity is that II rise to some value and 
thereafter remain constant and less than 2.6. 

Equation (17) describing the development of II shows that 
if //<(l.4+ (— ZV/A), then dH/dx^> 0 and H tends toward 
1.4+ (— Dr/A). If Z7>1.4+(— D t /A), then equation (17) 
shows that dll/dr<i 0 and therefore II always tends toward 
the value 1.4+ (— D T /A). 

If constant values of D r are considered as possible con- 
ditions for determining the velocity distribution, then the 
momentum equation shows that for a given velocity change 
D t should be as large a negative number as possible in 
order to reduce the surface for skin friction. (The rapid diff u- 
sion results in airfoils with relatively longer regions of constant 
maximum velocity and therefore larger airfoil circulations.) 
For reliability of the boundary-layer equations, it is desir- 
able to maintain P r > — 0.01. Therefore a possible condition 
is D t =— 0.01, the largest rate of diffusion with/=0.006534. 
The limiting value of H for this case is 2.14, a value that is 
probably, but not certainly, safe. Two values were therefore 
used in the velocity computations: D T = — 0.01, and D T = 
—0.006, which corresponds to a limiting value of II— 1.844. 

A third condition for prescribing the diffusion rate was 
also employed. According to reference 4, separation will 
result if 

D s = ~ R-°™=^ 7 ~ ™ (18) 


attains the value of —0.0013. The third type of velocity 
distribution computed was for the maximum diffusion rate 
with both D t > — 0.01 and Z? s > —0.00115; D s was arbi- 
trarily assigned with the condition that it bn greater than 
the separation value —0.0013. 

Velocity distributions in region of diffusion. — For the 
cases where D T is constant, equations (12) and (17) are used 
to eliminate the differential of x, and to obtain 


_Dt dll 

U~ A e H-IA+D t /A 


(19) 


where 


0-K-^)] 


( 20 ) 



is the exponential integral function (tabulated in reference 8), 
having the principal part of the integral for y ^> 0 (that is, 
for H<IAA—DtIA). The constant of integration for U— U 0 
was chosen at H=H 0 — 1 .4. This value of II and U 0 will not 
be included on any point of the airfoil if laminar diffusion 
precedes turbulent diffusion. 

Equations (12) and (13) give 


f-S [i-(T+ 14 - w)] T7 -l (»-»-*+$) 


dU 

U 

( 21 ) 


For the factor dU/U of the second term on the right side, a 
substitution is made from equation (19) to obtain an equa- 
tion that may be integrated to 



The arc length is then given by 

x-x 0 _ r z dx = ± 

To J 1 o Ta Dt] U 0 T q Tj 


(23) 


The curves for the boundary-layer parameters in the tur- 
bulent-diffusion region are shown in figure 4 for D r = —0.010, 
D r = — 0.006, and maximum diffusion rate with D T > — 0.01 
and D s >- 0.001 15. 

In solving the boundary-layer equations (13) and (17) 
(with /=0.006534) with the condition D s > —0.00115, the 
value of D t rises with R. The limitation D r >— 0.01 was 
added at the high Reynolds numbers in order to keep D T in 
the range for which the friction law /=0.006534 is valid. 
The solution of the equations in the region of variable D T 
was obtained by approximating each portion of the curve 
with the solution for Z? r =constant=mean value of D T in 
that step. 

RELATIONS BETWEEN VELOCITY DISTRIBUTION AND BLADE REYNOLDS 

NUMBER 

The relations between the velocity distribution and blade 
Reynolds number U r Xjv are developed. Some of the results 
are given in graphical form, which permits the use of the 
generalized velocity-distribution curves developed for the 
laminar and turbulent boundary layers (figs. 2 and 4). These 
curves will give the velocity distribution in the region of 
diffusion if the scale constants L h U a , and T 0 are known; 
the remainder of the surface is determined by the constant 
x t or x t . 


which integrates to 
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fc) Maximum diffusion rate Dv>— 0.01 and Ds>— 0.001 15. 

Figure 4— Concluded. Turbulent boundary-layer grosrth. 

In order to assist in finding the appropriate constants, the computing fig. o is described in appendix B.) The curves 

quantity x t fX (or x,fX for completely turbulent diffusion) is are to be applied in the following manner: 

plotted in figure 5 as a function of the Reynolds number The Reynolds number L'/Xjv and velocity ratio U r jU m at are 
UfXjv and the velocity ratio U f /U max - (The method of selected. Examination of figure 5 will indicate whether the 



Figure 5.— Variation of length ratios xrfX and rJX with blade Reynolds number. f?i,-250. 
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(b) Z>r« -0.006. 



(c) Maximum diffusion rate with Dr> — 0.010 and Z>$^— 0.00115. 

Figtjp.k 5.— Concluded. Variation of length ratios ujX and xJX with blade Reynolds number. Air=250. 


diffusion is laminar, mixed, or turbulent. If the diffusion 
is laminar, figure 5 gives the value o fz t a nd then £{=0.441 x h 
which establishes the scale for the use of figure 2. An- 
other procedure would be to determine {X-<r.i)/Li from 
figure 2 and solve for Zj, using £{=0.441 zj. 

If the diffusion is altogether turbulent, figure 5 determines 
x t and then .. . 



where (X—? 0 )/T a is given in figure 4 as a function of U f fU 0 
and Uq— Umax • 

The case of mixed diffusion is considered to bo of no tech- 
nical, significance in this application because of uncertainty 
in locating the position of transition. The effect of Reynolds 
number is therefore considered only in regions of purely 
laminar or purely turbulent diffusion, although x } /X and 
$/U f X were computed for Ri r = 250 and £ r = — 0.01 in the 
mixed-diffusion region. 
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Suction-surface potential. — A further step in computing 
airfoil performance is the integration of the velocity curves 
to obtain 



a dimensionless form of the potential on the suction surface. 
Results of the computations are shown in figure 6. The 
comparative magnitudes of the potentials for various values 
of D r (or D s ), UfjU mttx , and the Reynolds number are signif- 


icant because in the system used for estimating blade circu- 
lation, an increase in suction-surface potential is usually 
accompanied by an increase in blade circulation. 

The obtainable potential is highest with the greatest diffu- 
sion rate (compare figs. 6(a) and 6(b)) because the maximum 
velocity is maintained along the largest portion of the airfoE 
surface. The most pronounced effect of Reynolds number is 
the increase in potential shown in the change from laminar 
diffusion to turbulent diffusion (fig. 6(a)). The increase in 
potential begins at a Reynolds number in the vicinity of 
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(c) Maximum diffusion rate with Dr^— O.OlO and Ds>— 0.00116. 

Figure 6. — Concluded. Variation of suction^urface potential with blade Reynolds number. i?ir*250. 


6 8 /x/O 7 


80,000 to 90,000 and for C///t7 mo i>0.6 is fullyi realized at the 
blade Reynolds number approximately 260,000 for 
D s >— 0.00115 arid = — 0.006 and approximately 180,000 
for D t —-~ 0.010. This change increases with decreasing 
diffusion- velocity ratio U r /U max . Although there is no effect 
of Reynolds number on suction-surface potential for diffusion 
with a completely laminar boundary layer, the effect of 
increase in Reynolds number for turbulent diffusion is a 
decrease in potential for c onstan t values of D T . Am insignif- 
icant effect is shown, however, when the criterion of Kalikh- 
man for separation is applied (fig. 6(c)). Because at large 
Reynolds numbers the condition D s — — 0.00115 is replaced 
by D t = —0.01 , the curves so constructed approach the curves 
for D t =— 0.010 in that region. (The large Reynolds num- 
ber condition is equivalent to the assumption that the extent 
of the laminar boundary layer approaches zero.) 

The case of infinite Reynolds number is of particular signifi- 
cance because it is the most conservative design with 
turbulent boundary layer for constant diffusion coefficient P r . 
The potential at large Reynolds numbers was therefore 
computed and is shown in figure 7 for values of D T = — Q’.OlO 
and —0.006. The curve for Z? r = — 0.010 at high Reynolds 
numbers is practically the same as that of D s > —0.00115 for 
all Reynolds numbers because the potential for the D s con- 
dition shows so little variation with Reynolds numbers. 

The trend of curves shows that highest potential occurs at 
a Uf/Umai of approximately 0.58 with very little variation 
between 0.53 and 0.63. . Therefore, if a design is considered 
with a preassigned value for the velocity diagram and Uj, 
negligible gain occurs in using higher peak velocities than 
Umax** 1.6 Up For diffusion with a laminar boundary layer, 
figure 3 shows negligible increase unobtainable potential for 
peak velocities U max higher than 1.25 U f . 

Deviations from design conditions. — A knowledge of the 



Figure 7.— Relation between suction-surface potential and velocity ratio for blade Reynold* 
number approaching infinity. 


range of operation of any cascade of airfoils is essential in 
determining its suitability for application. For the sets of 
blades considered herein, the effect of change of Reynolds 
number is nonexistent with purely laminar diffusion. A 
very slight and variable effect is shown for turbulent diffusion 
with Ds> —0.001 15. If the limitations imposed by constant 
Dr are closer to the physical facts, less lift is obtainable at 
highetJReynolds numbers. Consequently, increase in Rey- 
nolds .number of operation would involve the danger of 
separation. A safe estimate of lift would therefore be made 
on the basis of. the potential for the maximum Reynolds 
number of operation, provided the minimum Reynolds num- 
ber of operation does not extend into the regions requiring 
laminar diffusion rates. If the maximum Reynolds number 
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of operation is taken to be three times the minimum, a very 
slight effect of design Reynolds number on potential is 
indicated even for the diffusion rates D T =— 0.010 and 
—0.006. The Reynolds numbers at which some diffusion 
with laminar boundary layer occurs depend on the value 
selected for R tT (in this case 250), on the diffusion” rate, and 
on the velocity ratio U^/ U f . Estimates of lift made on the 
basis of D t = — 0.010 and RN — would be in agreement 
with estimates based on Z? s > —0.001 15. If a blade designed 
on such a basis were used at low enough Reynolds numbers, 
the laminar boundary layer existing in the region of diffusion 
would separate because of too steep a pressure gradient. 

If a rough finish or a higher angle of attack than the design 
value occur in application of the blade, an increase in the 
value of T, over the design value may be expected. Then 

the initial value of will increase in the same 

proportion at the beginning of diffusion. Sample compu- 
tations show, however, that for velocity distributions based 
on constant D T , this initial increased value for D T is not 
maintained but drops to some value intermediate between 
the initial value and the value for which the velocity distri- 
bution was prescribed. If the value for separation of the 
form factor II is assumed to be 2.2, it would correspond to the 
limit for a constant value of D T given by 

D t = -0.0135 (7Z— 1.4)=— 0.0108 

If the design velocity distribution corresponded to 
D t =~ 0.006 (limiting 77=1.85), and then in operation T 
were increased by 80 percent to give an initial value of 
D t —— 0.0 10S, the resulting limiting value of D r would be 
greater than —0.0108 (less in absolute value) and the limiting 
value of II would be less than 2.2; thus, according to this 
criterion, separation would be avoided. If the design 
velocity distribution corresponded to D r —— 0.010, however, 
an increase of 80 percent in T would result in separation 
according to this criterion. 

ESTIMATE OF BLADE CIRCULATION AND SPACING 

EFFECT OF THICKNESS 

An approximate calculation is made for the blade circula- 
tion from the upper-surface velocity distribution by a modi- 
fication of the method for thin airfoils. The basic idea is 
that the airfoil-surface velocity is the sum of two components, 
one of which results from the basic thickness form un- 
cambered, and the second of which results from the curvature 
of the camber line. The suction-surface velocity of the thin, 
slightly' cambered, isolated airfoil is then given by 

U'=U'+AU'=V> (f5+ A T S-) 

where 

U' suction-surface velocity of airfoil 

U' velocity component result ing from uncambered thickness 
distribution 

A 77' velocity increment resulting from curvature of camber 
line 

V' free-stream velocity 

056646—51 10 


The ratios U'[V' and AU'/V' are independent of the mag- 
nitude of V f . A similar equation holds for the pressure 
surface, 

U b '=U , -AU'=2U'-U' 

where 77/ is the velocity on the pressure surface. If the 
suction-surface velocity is known, and if the velocity com- 
ponent U' resulting from the thickness can be determined, 
the pressure-surface velocity can be computed. 

In applying these ideas to the cascade, the principal change 
is the use of the distorted flow field in which the airfoil is 
immersed; the quantity V' is now regarded as a variable 
quantity. Near the leading edge, the velocity Tu is substi- 
tuted for V giving 

U t =v i(fsT) (26) 

where the _subscript i indicates the value near the leading 
edge and U/l V is the velocity rise over the free-stream 
value for the isolated airfoil. The value for U'/V' is a func- 
tion of the assumed thickness distribution. Similarly, near 
the trailing edge the free-stream value used is V 2 giving 

U,= V t (¥£) (27) 

At the point of the airfoil where the air has been turned 
halfway so that the flow direction is given by the angle a m 
where 

tan (tan a, -{-tan or 2 ) (28) 

where 

a angle between velocity vector and normal to cascade axis 
Subscripts: 

1 value at position upstream of cascade 

2 value at- position downstream of cascade 

the same procedure might be used as for the leading and 
trailing sections of the airfoil with the substitution of V m 
for V' in the equation for U m . The value of V m is computed 
from the continuity equation 

V L cos ai = V 2 cos a 2 =V m cos a„ t (29) 

which does not include the effect of airfoil thickness in block- 
ing the flow area and increasing the flow velocity. If the 
spacing s is not too small, the correction for the distorted 
flow that is applied is 



This formula is derived in appendix C. Then, 
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When the blade spacing is very small, however, it becomes 
more accurate to regard the flow between the blades as flow 
in a channel. From this standpoint, the flow velocity is 
increased over V m by the ratio of the flow area without blade 
thickness to the flow area with blade thickness. That is, 
the effect of thickness is given by 

where t m is the blade thiclmess.at the point x m where the air 
flow has been turned to the direction a m . This equation is 
applied only if the point of application of the lift force x m 
is within a range of values specified by 


^<x m <X- 


(32) 


This condition was selected because inside this range the 
effect of variations of inflow _and outflow direction on the 
velocity will probably be small enough to be neglected for 
the purpose of the calculation.. Thus, the flow is expected 
to be essentially that inside a long channel. 

Relations have been developed for computing TJ at three 
points on the airfoil. The family of airfoil cascades being 
considered will be restricted by assuming a second-degree 
parabola connecting the points (0, t7 4 ), (x m ,U m ), and (X,U r ). 
That is, _ 

vrrf+ Ba i+ B “(£)' -m 


where 


B «>• 




UJU, 


wr Bm 

UJU, 


Ui 

U f 


VflU, 


xJX (a :JX) (1— a :JX ) T (1 -xJX) 


(34) 


ESTIMATE OF CIRCULATION AND SPACING 

When TJ U U m , and Uf are assumed to be known, the curve 
for TJ can be determined when the value of x n is found. 
Because it is the center of pressure, x m may be found as the 
weighted mean value of a:; the weight assigned to each value 
is the local blade loading. The definition of x m may be 
reduced to the form 


J] xd{<p~ <p p ) 


(35) 

where • 

<P potential on upper surface measured from leading edge to 
a point a: ... ■ ■ 

r airfoil circulation, ($— $ p ) 

$ value of <p at trailing edge ~ 

The quantity, 

may be found from figures 3, 6, or 7. The quantity 

j xd(/>— S xUdx=UfX 2 m ........ 

Jo Jo 


may be found from figure 8 for the conditions of laminar 
boundary layer and very high Reynolds number. The 
curves for D r =— 0.010 are approximately correct for 
As>— 0.00115. ...... _ 



For calculation of the potential <p v> the high camber, of 
airfoils in.caseade must be considered because it results in a 
substantially different length of the suction and pressure 
surfaces.unlike isolated airfoils, which are nearly straight. 
The effect of the difference in length is approximated by 
setting ' 


(36) 


x p =$x 
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where /S is some constant nearly equal to unity. The circula- 
tion is thus 


r=J*Udx-j o Xr Up dxp= j o X [U-p( 2 U-U)\dx 
F—§—( 2 [}^ X Udz—p$ S j=$—$p (37) 


where 


$ p =2 pj* U dx—p$ (37a) 

When equation (33) is used for U, there are obtained 


curve for U (equation (33)) determined. Furthermore, 
equations (37a) and (38) determine the potential so that 
the dimensionless circulation TjUj-X is then 


r 

u f x~u f x u f x 


(44) 


For potential flow, 

r=s(F l sin at — Fj sin a . 2 ) 

The velocity Fi may be eliminated by the continuity equa- 
tion (24) resulting in 


and 

with 


f% d<p p =pj X x Up dx=t 3 ^ X x (2 U- U) dx 

l Z V zdt-U'X* (i §4-^4*") (39) 


The factor /3 is estimated on the assumption that the airfoil 
suction surface and pressure surface are circular arcs inter- 
secting at the leading and trailing edges. The angle of 
turning X of the mean line is assumed equal to the air 
deflection: 

X=a 1 — aj 


X V 3 cos ai (tan au — tan a 2 ) fA! .\ 

l~W f r C45) 

u f x 

which determines the blade spacing. 

COMPARISON WITH EXACT PROCEDURE 

For the purpose of judging the accuracy of the method of 
estimating the pressure-surface-velocity distribution, a com- 
parison was made with velocities computed for potential flow 
by an exact procedure. The data were obtained from refer- 
ence 9. The comparison is shown in figures 9(a) and 9(b). 
In both cases, values of U'/V' were selected to fit approxi- 
mately the NACA 65-series isolated airfoil data (reference 10), 
which gave 


The turnin g of the suction surface X is then computed from 

(40a) 

Then the turning of the pressure surface X„ is computed from 


• t X cos (X/4) 

SU1 4 X 4 sin (X/4) 


tan^-tan^— 2 
The ratio of lengths is 

Xp = \p sin (X/2) 
p X X sin (X„/2) 


(40b) 


(41) 


when the values for the integrals (38) and (39) are substituted 
in equation (35), 


C< 31 ( ^) 3 + C« + C m (j) + T2T7' = 0 (42) 


where 


10(3) = 

_u c +u ; 

■ fl+t S\ $ 


L = 

- 2 Ur 

V 2/3 ) UfX 


n+p 



5 U t 

\ 2/3 y 

1 UfX 1 \ 

. 2/3 ) m 3 U f 

"6 U r 


C m = 


_ 5 U t j 1 Uf /' 1+fl ^ m —— 


12 U f ~l U r V 2/3 J m 6 U f J 


(43) 


With the values for Ut, U m , and U t equation (42) may be 
solved for x m , equation (34) for B M and B i2) , and finally the 


fft ") 

^7=1.025 

ff f t 

^-=1 + 1 . 35 ^ - 
77 ' 

y?r= 0.950 


(46) 


By using equations (26), (27), (29), and (30) with these 
values, U f /V m , Uf/V m , and UJV m were computed. Next /3 
was obtained by use of equations (40) and (41). Evaluation 
of the data for suction-surface-velocity distribution yielded 
3>/V m X and m(Uf/V m ). It was now possible to substitute all 
these values into equations (43) and solve for x m fX by use of 
equation (42). The solution was then continued in the 
manner indicated in the previous section. 

Pressure-surface velocity was then plotted using the 
equation 

rj u+u, 

U 2 


The shapes of the estimated and exact velocity distribution 
on the pressure surface are quite similar. The solidity com- 
puted from the estimated velocities is, however, 6 percent 
higher for figure 9(a) and 4 percent lower for figure 9(b) than 
the exact value. These values then give an indication of the 
order of inaccuracy of the method of estimating solidity. 
The error in pressure-surface potential was 6 percent for 
the airfoil of figure 9(a) and 2 percent for the airfoil of figure 
9(b). 
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O 20 40 , 60 80 IOO 

Chord, percent 


(a) NACA 85- (12)10; c/i-1.00; w -83.0”; r/V»c=0.735; <r„-63.4°; estimated c/r-1.06. 

(b) NACA 65— (8)10; c/»=I.OO; at=3.7°; r/V»e-0.305; a B — 0.8°; estimated c/s-0.90. 

Figure 9.— Comparison of estimated and oomputed pressuro-surfaeo- velocity distribution. 

EFFECT OF SEVERAL PARAMETERS ON REQUIRED CASCADE 

SOLIDITY 

The relations developed in the previous section are applied 
to study the effect of several design and operating parameters 
on required cascade solidity. 

COMPUTATIONS 

The data given in the problems solved consisted of: 

(1) Reynolds number UfX/v (either very high or very low) ; 

(2) diffusion coefficient D T (either —0.010 or —0.006 for a tur- 
bulent boundary layer) ; (3) trailing-edge loading ( U f — U/)/V 2 
(assumed to be 0.1 because cascade data indicated this value 
as attainable); (4) velocity ratio U max /V 2 ; (5) velocity- 
diagram components Vi, V s , a u and a 2 from which V m and 
a m are computed; and (6) velocity ratios Ut[V r , U m '/V', and 


U/ fV', which are determined from the assumed basic thick- 
ness distribution. In all cases the assumed values were 
U t '{V'= 1.025 and Uf/V —0.95, but U m 'jV' was varied to 
obtain various values for the thickness. All the assumed 
quantities were varied in the examples to show the effect on 
the required cascade solidity, the velocity distribution, and 
the thickness. The sequence of computation is: 

(1) Compute 


Y~ T77— 1.°25 


0.95 


y~Y T ~ 

1 05 

Vi vEKrv, y 


U f 


Ei 

v 2 


u m 


U K 


v 2 


(2) From Uj/U maz , the assumption of either laminar or 

turbulent diffusion, and the desired rate of diffusion, figures 2 
to 8 are used to find the suction-surfacc-velocity distribu- 
tion and the values for $/U,X and m. .ZT 

(3) Values for the solidity X/s and the pressure center 
x m /X are assumed. The thickness is then computed from 
equations (30) and (46) if either z m /s<0.5 or r m /s>JV/s— 0.5. 

If 0.5<r m /s<Z/s— 0.5, then equation (31) is used. 
The factor |8 is then obtained from equations (40) and (41). 

(4) The quantities C°\ (7 t2 \ (?“>, x n , B ll) , B™, $/U f X, I 
T/UfX, and X/s are next computed from equations (43), (42J, 
(34), (38), (44), and (45). 

(5) Steps (3) and (4) must be repeated if the solidity X/s 
varies considerably from the assumed value or if, because of 
an incorrect assumption for x m , an incorrect choice was made 
in the equation for the computation of t m /X and /3. 

RESULTS 

The results of the computations are shown in the following 
table and in figure 10. The table gives the assumed values 
for the design parameters and the computed solidity; whereas 
figure 10 shows the details of the velocity distributions with 
sketches to indicate the cascade solidity and the airfoil 
thickness. The circles drawn on the camber lines indicate 
the location of x m and the airfoil thickness t m . 


Case 

Boundary layer 


“i 

(deg) 

<*I 

(deg) 

E 

v,-V, 

X 

X 

s 

V, 

a 

Turbulent 

-0. 010 

45 

38 

1.015 

0.10 

0.10 

0.30 

b 

Turbulent.- 

- .010 

45 

38 

1.312 

.10 

.10 

l&EB 

c 

Turbulent 

- .010 

45 

88 

1.3X2 

.10 

.00 


d 


- .010 

45 

38 

LS12 

.00 

.05 

.05 

0 

Turbulent 

- .010 

45 

-45 

1.615 

.10 

.20 

l.G 

f 


- ,010 

45 

-45 

1.615 

.10 

.16 

l.S 

s 

Turbulent 

- .010 

45 

-45 

1.812 

.10 

.u 

1.8 



— . 010 


-45 

1.312 

.10 

.18 

1.0 

i 


— .000 

45 

38 

1.615 

.10 

.10 

.42 

i 


- .098 

45 

38 

1.312 

.10 

.00 

.50 

k 

Laminar 

r- . 06818 

45 

38 

1. 815 

.10 

.08 

1.0 

1 

Laminar 

- ■. 00818 

48 

-48 

1.312 


.03 

3.1 

m 

Laminar... 

- .00018 

0 

-45 

1.312 

.10 

.09 

1.7 
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fa) Case (a) Dr—0.010; RA r -><»; «,=45°; ^=38°; (b) Case (b) Dt= -0.010; R.V-x=; <n=45°; ai=3S°; (e) Case (c) flr= -0.010; ai-45 0 : a ! =3S' > ; 

t<*!X= 0.10; X/i=0.31. tJX=0.W; J57»=0.64. t m IX= 0.0G; A?i=0.50. 

Cd) Case (d) Dr- -0.010; f!A r -»<=; <n=4S°; ai=3S°; (e) CaseCe) £>r=-0.m.O; ai-45°; m=-4S < ’; Cf) Case (0 Cr= -0.010; RA r -+<=; c«-4S 0 ; aj=.-4i g ; 

f«,/X=0.06; X/t=0.65. tjx=0.20; X/t=1.6. f „/X=0.1S; .X7i=1.3. 

Figure 10.— Examples of velocity distributions. 








Veiooify ratio Velocity ratio 
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fe) Case (g) -0.010; RN-»<»: 01-45°; as- -45°; 

(»/A«0.11; A/s=1.8. 

0) Case (J) Dr= -0.000; RlV->®; on-45 0 ; oi=38°; 
WA-0.00; A/«- 0.58. 


(h) Case ft>) iJr=-0.010j fiW4«; oi-0°; oj=-45°; ffi Caso (i) -Dr- -0.005; K,V->»; oi-46°; oj-38°; 

_<„/A-0.18; A/«=1.0. U'A-0.10; A/j-0.42. 

(k) Cose (k) Di— 0.00618; RN< 70,000; oi=45°; o a =-38°; 0) Cass 0) Dt— 0.00618; 1?A<70,000; oa-45°; ai--45°; 

(WA-0.08; A7«-1.0. ■ WA-0.08; A?s-3.I. 


Figuee 10.— Continued. Examples of velocity distributions. 


l 
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0 .2 .4 .6 , .8 1.0 

*PC \ 


\ 

(m) Case (m) Dl— 0.06618; KA r <70,000; oi=0«; <n— 45°; f«/X=0.l)9; Xfr-1.7. 

FlGl’RB 10.— Concluded. Examples of velocity distributions. 

The effect of variation of maximum surface velocity 
on solidity Xfs of a compressor-type blade is shown 
by comparison of cases (a) and (b), where all design para- 
meters are the same except U max fVi. The blade with a lower 
value of U m axl V 2 is obviously better suited for operation with 
a compressible fluid with a high Mach number, but has the 
disadvantage of requiring a higher solidity (twice as many 
blades). 

If in case (b) the thickness is reduced from 0.10 to 0.06, 
the cascade corresponding to case (c) results, having only 78 
percent of the solidity of case (b). This decrease results 
from the higher aerodynamic loading of the t hinn er blade 
because of the lower pressure-surface velocities caused by 
the thickness. 

In cases (a), (b), and (c), the trailing-edge loading corre- 
sponds to an air flow that does not satisfy the Kutta con- 
dition. The actual blade having the prescribed velocity 
distribution would have an extension from the point x/X=l.O 
to larger values in which the loading would drop to zero. 
Similarly, the initial parts of the velocity distributions do 
not include the leading edge; that is, the velocity distribu- 
tions shown do not extend over the entire airfoil, but do in- 
clude practically all the aerodynumieally loaded region. A 
velocity distribution was computed for zero trailing-edge 
loading (case (d)) in order to observe w'hat sacrifice would 
be involved in avoiding separation at the trailing edge and 
attempting an airfoil design that would satisfy the Kutta 
condition. Comparison with case (c) shows that the solidity 
is about 1.3 times as great. The effect would be less in the 
case of reaction blading where the low inlet velocity results 
m large loading in the initial part of the airfoil and a relatively 
smaller effect of trailing-edge loading. The effect of decrease 


in trailing-edge loading is also less for blades of small diffusion 
ratio Umaxl Uj-, because in such cases the suction-surface 
velocity is maintained at its maximum value to a point very 
close to the trailing edge and, as a consequence, the loading 
is also maintained over most of the airfoil. 

For an impulse blade (V a =T’i) of large turning (turbine 
type), the velocity distribution is very much the same as for 
a compressor blade with nearly the same velocity curves U 
and U, although the thickness will vary radically because of 
the low value for V m - For example, compare cases (a) and 
(e). 

If the turbine-blade thickness is reduced to the value closer 
to the compressor blade, the circulation is increased (case (f)), 

as might be expected. The large diffusion ratio on the 
pressure surface, which practically insures the existence of 
a separation bubble there, should also be noted. 

The effect of decreasing the maximum surface velocity U m 
is shown in case (g), which has lower blade circulation even 
though the thickness was decreased from 0.15 to 0.11. Solid- 
ity increased 39 percent over the value for case (f). 

For an inlet guide vane turning the air from an angle of 0° 
to an angle of — 45° (case (h)), the velocity distribution 
changes from the impulse cascade (case (g)) because of the 
large decrease in the inlet velocity. Even with the increased 
thickness, local velocities on the pressure surface do not in- 
volve high diffusion rates because of the over-all increase 
in velocity (Fi/Fi= 0.707). 

The effect of reducing the rate of diffusion from that of 
eases (a) and (c) (D T =— 0.010) to a lower rate {D T —— 0.006), 
while maintaining the same values of the other design 
parameters, was determined by computing the corresponding 
cases (i) and (j), respectively. The positions at which the 
UfV 2 values are equal for the two cases are where the values 
of a are the same. Some chord wise shift occurs in the loca- 
tion of a given value of a when changing from case (a) to. 
case (i) because of the shift in the loading distribution. 
Consequently, the location of x m shifts and the thickness 
distribution changes. The last half of the blade (case (i)) 
does very little turning so that the solidity increases from 
0.31 to 0.42. The usefulness of this part of the blade could 
be greatly increased if a larger turning of the air were in- 
volved because the U/V 2 curve would rise to higher values in 
the first half of the airfoil and drop to lower values in the last- 
half, thus shifting the load backwards on the blade surface. 
The change in velocity and solidity from case (c) to case (j) 
is slight because the rate of diffusion has very little effect on 
suction-surface potential when the diffusion is small. (See 
fig- 7.) 

An effect of blade Reynolds number U r X[v is shown by 
comparison of cases (a) and (k). In case (k), the expansion is 
entirely laminar (low Reynolds number) ; whereas in case (a) 
the extent of the laminar region is assumed zero (high 
Reynolds number). Even though the laminar-layer blade 
is only 6 percent thick, the solidity (1.6) required is five 
times as large as for the blade with a turbulent boundary 
layer. This loss in circulation results from the fact that 
laminar diffusion must be begun so much earlier when the 
boundary layer is thinner. On the pressure surface of the 
blade can be observed a region near the trailing edge where 
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the rate of diffusion increases' to a value that will certainly 
induce separation. A better aerodynamic design could be 
obtained with a thinner blade, which would be necessaiy if 
the turning were greater. Some improvement would also 
probably result with a less simple mean velocity curve. The . 
conclusion may safely be made that it is extremely difficult, 
if not impossible, to obtain a compressor blade of reasonable 
thickness and turning, if flow separation is to be avoided with 
a laminar boundary layer. 

The design of an impulse blade with laminar boundary 
layer appears to be entirely feasible (case (1)) as a consequence 
of the low value for TJm/W, which is a result of the.low value 
for VJV,. A. reaction blade with laminar boundary layer 
also is feasible. (See case (hi).) The. percentage increase 
in solidity for cases (1) and (m) over that for cases (g) and 
(h) is much less than the percentage increase in solidity for 
case (k) over that for case (a). As a general rule, cascade 
solidities for laminar and turbulent diffusion can be expected 
to be more nearly alike when diffusion is small (UmazlU/— >1.0) 
than when diffusion is large (Z7 mM /Z7/-»2.0) because of the 
smaller difference in suction-surface potential. (See fig. 6 (a) .) 

COMPARISON OE EXPERIMENTAL WITH INDICATED 
ATTAINABLE SOLIDITY VALUES 

In order to estimate the degree with which blades in use 
approach the maximum attainable circulation as indicated 
by the present method of analysis, cascade solidities were 
computed for conditions comparable to those for which 
performance data were available. It was desired to compare 
the solidities of the tested blades and those corresponding to 
the proposed velocity distributions on the basis of operation 
at veiy low Mach numbers, very high Reynolds number, the 
same values for a 3 , maximum thickness t/X, Umaz/Vt, and 
UfjV-i. Data used were obtained from. references 11 and 12. 
In none of the comparisons shown is it possible to evaluate- 
cascade efficiency, although for every case the comparison 
was made on the basis of the optimum operating condition of 
the cascade as stated herein. The comparisons are shown in 
the following table: " " 


Ref- 

erence 

Blade 

designa- 

tion 

(deg) 

at 

(deg) . 


Experi- 

mental 

xa 

.. . Estimated X!» 

Dr*— 0.010 

Dr=— 0.006 

11 

65- (IS) 10 

45 

16 

1.699 

1. 003 

*1.08 


11 

65-02)10 

45 

23.8 

1.602 

1.034 

1. IT 

1. 60 

11 

65- (4) 10 

45 

87.0 

1. 30CT 

1.013 

.. .78 

.88 

11 

65- (0)10 

45 

42.9 

1.186 

1.003 

.36 

.38 

12 

64-(A)QG 

0 ‘ 

-45.8 

1.180 

1. 538 

.. 1,12 


12 

64-(B)O0 

0 

-62.0 

1* 133 

1.534 

1.20 



0.0664. 


All of the cases are. directly comparable as previously 
described except the 65-(18)10 blade, for which it was im- 
possible to find a solution with 10-percent thickness, indi- 
cating a strong probability of separated flow on the blade. 
A solidity was computed with, the thickness reduced to 6.64 
percent. 

The 65-(12)10 blade has a’Tower solidity and therefore a 
higher circulation than either of the estimates for a 
maximum-circulation blade. Possible causes for this discrep- 


ancy aref (l) slight flow separation from the surface of the 
65-(12)10 blade, and (2) the assumption of a parabolic 
curve forJ7, which may involve some sacrifice in attainable 
circulation. For lower diffusion (65- (4) 10 and 65-(0)10) 
and for reaction blades (64-(A)06 and 64-(B)06), the esti- 
mates indicate that higher circulations are attainable. The 
64-(A)06 and 64-(B)06 blades have such high values for 
Uf/Umaz that there is no substantial difference when the lower 
diffusion rate D T =— 0.006 is used in the solidity estimate. 
It seems reasonable to expect that higher blade circulations 
are obtainable with blades, other, than the 65- and 64-series 
for low pressure rise and for pressure drop. The 65- (0)10 
blade is particularly subject to improvement. 

- CONCLUSIONS 

The analysis of the limitations on the circulation about 
blades in cascade indicates that under the assumption of the 
several criterions for separation and the avoidance of local 
separation of the flow, the following conclusions can be 
drawn: 

1. For a preassigned maximum velocity and suction- 
surface trailing-edge velocity, the suction-surface potential 
and airfoil circulation increases with the increase in permis- 
sible diffusion rate. 

2. For_a laminar boundary layer, Loitsianskii’s equation 
determines blade circulations independent of the Reynolds 
number. For a turbulent boundary layer, Kalikhman’s 
criterion for safe diffusion (D s >— 0.00115) indicates very 
slight changes in blade circulation with variation in Reynolds 
number. If a constant diffusion rate (D r = — 0.010 or 
— 0.006) is used, however, some slight decrease in attainable 
potential with increase in Reynolds number is indicated. 
A very marked increase in obtainable lifts occurs when the 
design is changed from that for a laminar boundary layer to 
that for ..a turbulent boundary layer. This effect decreases 
with decrease in the diffusion ratio U mttX IUf or U ma t/Vi. 
With transition at momentum-thickness Reynolds number 
f?— 250, the change begins to occur at a blade Reynolds 
number RN between 80,000 and 90,000. For completely 
turbulent diffusion with U f jU ma X>0.6, the blade Reynolds 
numbers are 260,000 . (Kalikhman’s criterion and 
D r =-0;t)06) and 180,000 (D r =-0.010). 

3. Suction-surface potential and consequently the airfoil 
circulation increases with increasing maximum velocity 
ratio UmazlVi or U max jU t ; however, negligible gain in cir- 
culation is obtained with U max !Uf greater than 1.6 with tur- 
bulent diffusion and greater than 1.25 with laminar diffusion. 

4. An airfoil designed for laminar boundary layer may be 
operatecLat. any Reynolds number. An airfoil designed for a 
turbulent boundary layer has a lower Reynolds number 
limit for operation without separation. This limit is a func- 
tion of the velocity ratio U max IU/, the diffusion rate, and the 

■ conditioner transition. For R = 250 and diffusion velocity 
ratios l^u/UfK 1.67, the lower limit of RN is equal to or 
less than 260,000 for small diffusion rates (D T =— 0.006, 
D s >— 0.00115) and equal to or less than 180,000 for larger 
diffusion rates (D T —— 0.010). 
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5. Improvement in airfoil circulation always results from 
adding loading at the trailing edge, but the improvement is 
less with decrease in over-all diffusion (decrease in VJVt) 
and with decrease in surface diffusion (decrease in U max ! U f ) . 

6. Change in suction-surface-velocity distribution result- 
ing from change in diffusion coefficient or character of the 
boundary layer has a large effect on blade circulation except 
when the diffusion is small. This small effect is to be expected 
in blades designed for high-speed operation and for reaction 
blades. 

7. For a fixed suction-surface-veloeity distribution, ob- 
tainable blade circulation decreases with increasing blade 
thickness. 


8. The large discrepancy between the circulation of the 
65— (18) 10 blade with the recommended velocity diagram 
and the estimated attainable circulation implies a strong 
possibility of separation of the flow from the 65-(18)10 
blade. For the recommended conditions of operation of the 
65-(4)10, 65— (0)10, 64-(A)06, and 64-(B)06 blades, esti- 
mates indicate that higher blade circulations are attainable. 


Lewis Flight Propulsion Laboratory, 

National Advisory Committee for Aeronautics, 
Cleveland, Ohio, September 1, 1949. 



APPENDIX A 


SYMBOLS 


The following symbols are used in this report: 
A 0.0135 


B*\B m 

C u \C a \C a) 

c 

D l 


D s 


coefficients of (x/X) and (x[X)* terms in equation 

for tl _ . 

coefficients of (x m /X), (x m jX ) 2 , and (x m jXf 
terms in equation for x m 
chord 

diffusion coefficient with laminar boundary 


, /LdU\ 

\vtx) 

(according to Kalilthm&n, separa- 
tion occurs when D s = —0.0013) 


D r 


Ei(y) 

F 

f 

H 

L 

m 

R 

RN 

s 

T 

t 

U 

U 

V 


diffusion coefficient with turbulent boundary 

la * er ’ (uif) 

exponential integral, 

parameter in approximate equations for laminar 
boundary layer 

coefficient of. friction, R i/& ^j 

form factor, {&*!&) 

generalized momentum thickness of laminar 
boundary layer, ( dR ) 

UJT 2 fo xUdx 

momentum- thickness Reynolds number, (Ud/v) 
blade Reynolds number, ( U f Xjv ) 
spacing of airfoils along cascade axis 
generalized momentum thickness of turbulent 
boundary layer, (OR 116 ) 
thickness of airfoil 

suction-surface velocity just outside boundary 
layer 

arithmetic average of suction- and pressure- 
surface velocities 

stream velocity ' ' 



X 

x 

a 


& 

r 

5* 

6 

X 

X 

v 

P 

'r 

$ 


<p 


V m corrected for thickness effect of airfoils in 
cascade 

total suction-surface length 
„ . arc. length on airfoil suction surface 
- angle between velocity vector and normal to 
: cascade axis 

ratio of pressure-surface length to suction- 
surface length, {xp/x) 
circulation around one airfoil 
displacement thickness of boundary layer 
momentum thickness of boundary layer 
_ turning angle of upper surface 
, turning angle of camber line 
kinematic viscosity of gas 
* gas density 
37 surface shear stress 

complete velocity p 

velocity potential, 
leading edge 


lotential, (j> & ) 

(!>*) 


measured from 


Subscripts: 

/ 

i 

l 

m 

max 

p 

t 

tr 

0 . . \ .. 


final value near trailing edge 
initial value near leading edge 
value at beginning of laminar diffusion process 
value corresponding . to mean of upstream and 
downstream velocity vectors 
maximum value 
value on pressure surface 

value at beginning of turbulent diffusion process 
value at transition 

value corresponding to 11= 1.400 in turbulent 
boundary layer 


1 upstream value 1 

2 downstream value 


Superscript: 


pertaining to isolated airfoil 
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APPENDIX B 

METHOD OF COMPUTING FIGURE 5 


For diffusion with a laminar boundary layer, equation (8) 
determines for any given velocity ratio U f jU ma x (equal to 
Uf/Ui) and the blade length X the value for x t . Then 

Z{=0.441xt 


For the case of diffusion with a partly laminar and partly 
turbulent boundary layer, a number of values were selected 
for 



and 

T tr =6 tr {R tr y' 


The initial condition for H is 


Htr — 1.4 


Dt.IT 

0.0135 


(reference 3), where 


the Reynolds number at the beginning of laminar diffusion. 
The units of velocity and length may be taken as U t and x u 
respectively JjfcAlso 

Li=0A4lX[ 

which establishes the scale of the laminar-diffusion region. 

8 U 

At the transition^point, R tT =— — -=250, so that for L tr 
there is obtained from equation (6d) 



Ur 


and 

Uct 

IDl 

. U R *Y +tDk 

Timm \Rl) 


i rr 

Because 


s 

II 

then 


„ Ut 



&tr — ‘ p 
It It 


For this value of H tr , figure 4 identifies U tt JUo and Rit/Ro, 
thus establishing the scale for the velocity U and the basic 
length To- By proceeding to the final velocity U r , the length 
X—x tT is established, and RN= U r Xjv may be computed. 

If the diffusion is entirely turbulent, then by assuming a 
value for the Reynolds number at the beginning of turbulent 
diffusion R t , there are obtained (x tI and XJ tr are taken as the 
units of- length and velocity, respectively) 

»„-0.441 £ 

rp _ 0-44Ix tr 

1 ‘ r ~(R l r) SI * 

Dt.it — 0 

Htr=lA 

Dq D max 

T a =T t r 
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By assuming that the effect of the airfoil circulation, on 
the velocity at the point x m (where the flow has turned to the 
direction a m ) is taken into, account by the turning effect, the 
change in velocity from Iq to V m is given by 


APPENDIX C ' 

EFFECT OF AIRFOIL THICKNESS ON FLOW THROUGH CASCADES 

This approximate velocity distribution is seen from the 
form of the velocity equation to bo that of a doublet of 
strength 


V m =Vi cos'cti sec a m 

The effect of thickness is now separately evaluated as the 
approximate effect of a cascade of uncambered airfoils of 
stagger a m with no circulation. The thickness distribution 
is assumed similar to that of. a Joukowsky airfoil, which is 
discussed generally in reference 13. The potential function 
W for a uniform stream of velocity V flowing about a 
Joukowsky airfoil of small thickness can be shown to be 


W=V 


where 


(1 + 2e)c 2 /16 ' 

f+J- 


6 real constant determining airfoil thickness 

c 2 

f complex parameter defined by 2 =f+Ygji 

and the condition that f be a continuous, single-valued func- 


tion with lira 


coordinate in t 


£^-1=1, where z is the complex position 


» L ^ J 


he airfoil plane. The airfoil profile is given 


by the additional equation for a circle of radius (l+«) ^ 


and center—^: 


f+r=t 


When the central angle of the circle w has a value of.2/3 x, the 
airfoil has a maximum thickness ; at this point the velocity is 

U=V( l+2e) 

If the potential function is expanded in powers of 1/f, the 
individual terms may then in .turn be expanded in powers of 
l/a; the resultant approximation for large 2 is 


W 


= v(z 


jcV8_\ 


2 + c/8/ 


This expansion neglects terms in I/ 2 4 and higher powers. 
From this expansion, the velocity at large distances from the 
airfoil is - " 


dW 

dz 


—V- 


Vec 2 /8 

{z+cfS} 1 ' 


(U-V)c*/1Q 

(z+c/8) 2 


At a half chord above the airfoil 




the error in the 


(U—y\ c 2 

disturbance velocity j g c j ^ a is 26 percent; whereas the 

error in the entire velocity is (for f/=1.2 V) 0.8 percent. 
At a distance of 1 chord, the errors are, 10 and 0,1 percent, 

respectively. At 1 h chord, the values are 5 and 0.02 percent. 


Af=2x((7— F) jg 

placed in a uniform stream of velocity V. 

If a series of such doublets are placed on an axis with 
directionJe -< “’» and with locations 

z„±inse~ iam (n= 0, 1, 2, . . . ) 

a potential function IF will result with the equation 

,_!fr_U.f, ( -W + . - ■ ■ ^ )] 


w Mr 1 . y-, / r 

W ~ 2x \_z— Z b \2 — z„+in8e~ tam ^ z—z 0 —mser 


TTT Me iam , . 
W= cot | 


If the potential due to the central doublet at z 0 is removed to 
find the modification of the uniform flow in which the central 
blade is located, there results 

rJT M ( , rv(z— Z q )~1 ise-*"" ) 

~2ise~ iam ( C ° [_ ise~ tam J x( 2 — 2 „)) 

By expanding in powers of ( 2 — z 0 ) and by neglecting powers 
higher than 2, there is obtained the approximation 


IF= 


Mx 2 


fi2ia r 


~2tt 3 8 2 

with the complex velocity 

d\V = M x 2 
~3z 2 x 3? 


‘(z— Zo) 




When added to the mean flow, the component normal to the 
mean flow may be neglected, giving 

Tr T , . M x 2 „ 

^m,c — ' iiiT 2jr 

The doublet strength has been shown to be 

(U-V m , c ) 

If this equation is substituted in the equation for F m , c , and 
the equation 

V &' 

vzrv' - - 

used to eliminate V, . . 

Or, if V m ,c is eliminated, 

_ " (?' 

U m =V m , c 


“ 


u'lV' 

- 

1--I 
L 48 

( f ) 

>(?-0 

1 COS 2 a m 
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